clear all
load ParEst_eh0_time0_Mex_girls ParEst_eh0_time0_Mex_girls
load trans_time1_age0_eh0_comp1.out
Chute = ParEst_eh0_time0_Mex_girls; %Chute(12)=2000;
% Chute =[ParEst_eh0_time0_Mex(7), ParEst_eh0_time0_Mex(8), ParEst_eh0_time0_Mex(9), ParEst_eh0_time0_Mex(10), ParEst_eh0_time0_Mex(11), ParEst_eh0_time0_Mex(12),...
% ParEst_eh0_time0_Mex(17),ParEst_eh0_time0_Mex(18),ParEst_eh0_time0_Mex(19),ParEst_eh0_time0_Mex(20),ParEst_eh0_time0_Mex(21),ParEst_eh0_time0_Mex(22),...
% ParEst_eh0_time0_Mex(24)];
HandleLo = @Principal_A2_n_1;

D_mom=trans_time1_age0_eh0_comp1; % Read transition moments
Dij=D_mom';

Dnf_2	 =	Dij(9);
Dni_2	 =	Dij(10);
Dfn_2	 =	Dij(11);
Dff_2	 =	0;
Dfi_2	 =	Dij(12);
Din_2	 =	Dij(13);
Dii_2	 =	0;
Dif_2  =	Dij(14);

Lower =[Dfn_2,Din_2,Dnf_2,Dni_2,Dfi_2,Dif_2];
LoBo =zeros(1,13);
LoBo = [Lower, 0.1, 0.1, 0.1, 0.1, -10, -10, -10];
HiBo = [.81 .81 .81 .81 .81 .81 10 10 10 10 -4 5 0];
                 options = saoptimset( ...
                     'PlotFcn',{@saplotbestf,@saplottemperature,@saplotf,@saplotbestx,@saplotx},...
                     'TolFun',5000,'ReannealInterval',6, ...
                     'InitialTemperature', 10^2*[1 1 1 1 1 1 10 10 10 10 6 5 10], ...
                     'TemperatureFcn', @temperaturefast, 'StallIterLimit',1000); % 14* (default)
[x,fval,exitflag,output]  = simulannealbnd(HandleLo, Chute, LoBo, HiBo, options)
ParEst_eh0_time0_Mex_girls = x;
save ParEst_eh0_time0_Mex_girls ParEst_eh0_time0_Mex_girls
